Angular Dependences of Third Harmonic Generation from Microdroplets 



On 



O 



> 

O 

o 
o 

+-> 



c 

o 
o 



X 



J. KaspaiW"), B. Kramer J. P. Dewitz^i, S. Vajda^, P. RairouxW0, B. Vezin( Q \ V. Boutou^, T. Leisner^, 

W. Hiibner^, J. P. Wolf*°), L. Woste^, K. H. Bennemann( c ) 
(a) LASIM (UMR 5579), Universite Claude Bernard Lyon 1, 43 bd du 11 novembre, 69622 Villeurbanne Cedex France 

(b) Institut fur Experimentalphysik, Freie Universitdt Berlin, Arnimallee 14, 14195 Berlin, Germany 

(c) Institut fur Theoretische Physik, Freie Universitdt Berlin, Arnimallee 14, 14195 Berlin, Germany. 

(February 1, 2008) 

We present experimental and theoretical results for the angular dependence of third harmonic 
generation (THG) of water droplets in the micrometer range (size parameter 62 < ka < 248). 
The THG signal in p- and s-polarization obtained with ultrashort laser pulses is compared with a 
recently developed nonlinear extension of classical Mie theory including multipoles of order I < 250. 
Both theory and experiment yield over a wide range of size parameters remarkably stable intensity 
maxima close to the forward and backward direction at "magic angles". In contrast to linear Mie 
scattering, both are of comparable intensity. 
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Studies of optical microcavities have received increas- 
ing interest both from fundamental [jjj as well as applied 
research ^,3 . The quality factor of microcavities formed 
by liquid droplets can become as large as 10 6 . This leads 
to very high internal electromagnetic fields under reso- 
nant excitation. Accordingly, nonlinear optical effects 
including stimulated Raman scattering (SRS), third or- 
der sum frequency generation (TOSF) and in particular 
third harmonic generation (THG) in microparticles are 
enhanced by several orders of magnitude compared to 
flat surfaces and bulk media . 

These enhanced effects, in turn, make nonlinear op- 
tics a promising spectroscopic probe for the character- 
ization of microparticles. With respect to applications 
in this area a lot of interest is focused on determining 
size, shape, and refractive index from the measured in- 
tensities. With the help of the recently developed lasers 
generating ultra-short pulses of high intensities, higher 
harmonic generation can thus act as an additional source 
of information. 

Optical TOSF in particular THG in droplets was first 
reported by Chang and coworkers in many organic and 
inorganic liquids Hff). Theoretical models for TOSF and 
angle integrated SHG in microdroplets have been pro- 
posed by several authors || . The angular dependence of 
SHG and THG by particles with size parameter ka < 5 
(ka = 2ita/\ , a being the radius of the particle) has 
recently been investigated by Dewitz et al. using a non- 
linear extension of classical Mie theory || . 

Here we report the first observation of the angular de- 
pendence of THG generated by micrometer-sized water 
droplets ranging from 8 to 32 /im (62 < ka < 248), as 
well as new calculated angular distributions for droplets 
with radius up to about 10.5 /xm (ka ~ 81). Theory and 
experiment yield in fair agreement the following three 
remarkable features: 

-The THG scattering profile of droplets having a radius 



of a > 8 (jm is dominated by only a few intensity max- 
ima at magic angles close to the forward and backward 
directions. 

-The intensity and the angular distribution of the THG 
signal depends only weakly on the particle size and the 
central wavelength of the laser. 

-Theory and experiment find comparable intensity 
close to the forward and backward direction. 

Each of these findings is in strong contrast to linear Mie 
scattering and is of considerable relevance for the appli- 
cation of nonlinear optical spectroscopy to micro-particle 
sizing and detection. 

Our theoretical model for higher harmonic generation 
by microdroplets is similar to the formalism of linear 
Mie theory [|l0| . Therein the electrodynamic response of 
spherical particles is described by matching the incident, 
internal and scattered field with the help of the boundary 
conditions of their parallel components. The expansion 
coefficients of the scattered field can then be expressed 
in terms of the coefficients of the incident field. On 
this basis Ostling et al. || developed a model for angle- 
integrated higher harmonic generation by microparticles. 
This theory has been extended to angular dependences 
by three of the present authors Q . The main idea of the 
approach is that in a first step the incident field gener- 
ates a charge a = X)^Li where a^ n ' oscillates with 
frequency ncu. This picture is motivated by the anhar- 
monic oscillator model || . In a second step the nonlinear 
components cr*™) radiate the n-th harmonic field. In the 
central approximation of the theory o-(") is set equal to 
[o-W]™, where is equal to the radial component of 
the polarization known from linear Mie theory [^). 
The n-th power of er^ can easily be obtained by ex- 
panding a« in terms of spherical harmonics taking into 
account the appropriate selection rules for the angular 
momentum and keeping only terms varying with e 1TlLlJt . 
Using en") in combination with the boundary conditions 
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n- [Dq"j — Din ] — 47^(7( ,l ) for the radial components and 
n x [E^J — Ej"^ ] = for the parallel components and 
the appropriate expansion of the fields in vector spher- 
ical harmonics pj] j, one determines then the radiated 
intensities from the radial part of the Poynting vector 
n • (E x H)|. Here the subscripts "in" and "out" denote 
the fields inside and outside the droplet |l2j . 

In THG multipoles of order I < 3-ka+S (S w 20 accord- 
ing to the required precision) must be kept in the multi- 
pole expansions in order to get a fully converged series ex- 
pansion for the fields. Although conceptionally straight- 
forward it is numerically very time-consuming to obtain 
higher-order multipole contributions. Consequently our 
previous calculations were restricted to size parameters 
ka < 5. For small spherical droplets in this size range the 
theoretical treatment showed for SHG and THG an in- 
creasing number of maxima which where distributed over 
the whole angular profile. They result from an increased 
number of multipoles. The angles of the intensity max- 
ima exhibited strong variations upon small changes of the 
droplet size parameter ka. Moreover, the enhancement 
of the intensities scattered in the forward direction -well 
known from linear Mie scattering |l3| , |l4| ]-, was amplified 
in SHG and THG dominating the angular profile. In 
the present paper we focus on larger particles (ka > 62) 
obtaining quite distinct results. 

The experimental setup is given in Fig. [l]. Monodis- 
perse water droplets were generated at a repetition rate 
of 1 kHz by a piezo-driven nozzle. Depending on the 
operation condition of the nozzle, droplets in the size 
range between 8 /im and 32 /im could be produced. The 
droplet generation was synchronized to a Ti:Sa femtosec- 
ond laser system, which produces ultrashort light pulses 
(80 fs) with pulse energies up to 400 nJ at a central wave- 
length of 820 nm. The scattered light was collected with 
a quartz fiber mounted onto a stepper motor driven go- 
niometer, then spectrally dispersed into a UV spectrom- 
eter and detected by a photomultiplier tube. Polarizing 
beamsplitters could be inserted in front of the fiber to 
allow polarisation sensitive measurements. The setup al- 
lowed to analyze the radiation emitted from the droplet 
in an angular range between 18° < 8 < 160°. The angu- 
lar resolution was 1°. The shot to shot variation of the 
particle size was controlled independently by the analy- 
sis of the Mie scattering patterns obtained with a He:Ne 
laser and by direct observation of the droplets under a 
light microscope; it was found to be below 2 %. The sta- 
bility of the droplet generation and the temporal and spa- 
tial overlap of the Ti:Sa laser beam and droplet stream 
was continuously monitored by observing the scattered 
laser light with a microscope which was equipped with a 
CCD camera readout. At higher fluences of the Ti:Sa 
laser, the droplets disintegrated into a cloud of small 
droplets with a diameter of about 1 jim.. 

Since the optical response of the droplets to the ultra- 



short laserpulses is much faster than their fragmentation, 
the light scattering signal was nevertheless detected. It 
was checked experimentally that background signal from 
remaining small aerosol droplets could be neglected com- 
pared to the initial droplet signal. In order to determine 
the dominant channel of the nonlinear optical response 
we performed spectroscopic measurements. We find as 
the most significant feature a large third harmonic peak 
at 273 nm. It showed a cubic dependence of the incident 
laser flucncc, as expected for a third order process. The 
SHG peak at 410 nm is rather weak. This is an indica- 
tion for the good spherical shape of the droplets for which 
SHG is forbidden by symmetry. The THG light was s- 
or p-polarized as the incident light, the amount of the 
cross polarization was only 2%. Moreover, no stimlated 
Raman scattering was observed. This is comprehensible 
with respect to the short laser pulses used and demon- 
strates that THG is indeed the dominant nonlinear re- 
sponse channel. 

In Fig. H the experimental results for the angular de- 
pendence of THG of water droplets having a radius of 
a = 10.5 /xm (ka ~ 81) are shown for angles between 
18° and 160° in the ip = 0° plane. Two pronounced 
intensity maxima appear around 26° and 155°, whereas 
almost no intensity is detected at angles between 50° and 
130°. The result is in clear contrast to the angular de- 
pendence expected for linear Mie scattering from droplets 
with ka « 81. Instead of a very complex structure with 
several maxima of varying angle and magnitude only one 
broad maximum appears close to the forward and back- 
ward directions. 

Moreover we find that the angular dependence is re- 
markably stable upon variation of the droplet size: the 
scattering profile is reproduced for particles from a = 
8 /zm to 32 ^m. Only a small shift of the "magic angle" 
of maximum intensity occurs from 26° to 30° . Thus the 
angular dependence does not only show a strong reduc- 
tion of complexity but also a stabilization of the angles 
of intensity maxima over a wide range of droplet sizes, 
although the number of contributing multipoles varies 
strongly with ka. 

Most importantly we find that the magnitude of the 
intensity maxima close to the forward and backward di- 
rection is of the same order. This is in sharp contrast to 
linear Mie scattering where forward scattering is strongly 
enhanced. 

In Fig. H we show the theoretical results for the angu- 
lar dependence of the THG intensities of a water droplet 
with radius a — 10.5 (ka « 81) in the full 0-range 
from 0° to 180°, in the inset the corresponding polar plot 
is displayed. Dominating maxima appear close to the for- 
ward and backward direction at angles 8 w 10° (where 
no experimental data could be recorded) and 8 ss 158°. 
Those maxima show unexpected stability upon droplet 
size variation. The stability of intensity maxima in THG 
strongly contrasts with linear Mie scattering, where the- 
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ory yields a strong change of the scattering profile upon 
variation of ka in a small range around ka ~ 81 fT(| . Fur- 
thermore the magnitudes of the intensity maxima close to 
the forward and backward directions are of same order in 
THG as it is observed in experiments (see Fig. ||) . This is 
again in clear contrast to linear Mie scattering for which 
the calculated ratio of forward to backward scattering for 
water droplets with ka w 81 is around 6700 |l6[ . 

The stabilization of the angular scattering profile in 
THG as a function of ka results from dominating res- 
onances occurring within the laser bandwidth. These 
sharp resonance peaks are enhanced by six orders of mag- 
nitude as compared to the background radiation. The 
dramatic enhancement results from the nonlinearity in 
THG, which amplifies the field effects. Remarkably, our 
calculations yield that every resonance displays the same 
angular pattern. Experimentally these resonances are al- 
ways excited due to the finite laser bandwidth of ±5 nm 
resulting from the femtosecond pulses and the particle 
size fluctuation of ±2%. In our case these effects have 
been fully taken into account in the calculations. In the 
experimental fca-range (ka = 81 ± 2) four resonances are 
excited. 

In backward direction the position of the two peaks 
with the largest angles (156°, 160°) are found by both 
theory and experiment. For decreasing angle the agree- 
ment fades since the theory shows a more complicated 
structure and thus a stronger oscillation of the intensi- 
ties. The damping of the experimental structure results 
from the finite angular resolution which is of the order of 
1°. For intensities less than 130° both theory and exper- 
iment show no structure. 

In the forward direction theory and experiment agree 
in a less satisfactory way. Calculations predict a strong 
maximum at an angle of 10°. This, however, is not ob- 
servable with our existing equipment. Experimentally a 
maximum is observed at 26°. Work is in progress to in- 
vestigate the angular range also below 18°. 

In conclusion, the angular dependence of THG of wa- 
ter droplets with sizes from 8 /im up to 32 /im are studied 
experimentally and for sizes around 9 and 10 /im theo- 
retically. In good agreement both investigations show a 
reduced complexity of the scattering profile and a sta- 
bilization of the angles of maximum intensities close to 
the forward and backward direction. Moreover the mag- 
nitude of the intensities in the forward and backward 
direction are of the same order of magnitude. These re- 
sults are in remarkable contrast to the characteristics of 
linear Mie scattering in this size range. 

The angles of the maxima near the backward direc- 
tion show good quantitative agreement between theory 
and experiment. It was shown that the stabilization is 
an effect of the excitation of sharp resonance peaks with 
stable angular dependence independent of the size of the 
droplet. 

This stabilization can be considered in analogy to rain- 



bow scattering which occurs in linear optics for large 
droplet sizes and reflects the onset of a transition from 
Mie scattering to geometrical optics. The latter is valid 
for a plane surface resulting from the limit of droplets 
with infinite radius. Due to the incorporation of higher 
multipolcs in THG this limit is reached for smaller 
droplets than in linear optics. 

Thus THG opens up new possibilities to character- 
ize particles in the Mie-size range between Rayleigh- 
scattering and geometrical optics. Especially the strongly 
enhanced backward scattering yields new perspectives for 
the remote sensing of the atmosphere. 
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FIG. 1. Experimental setup 

FIG. 2. Measured spectra of THG intensity of a water 
droplet with a = 10.5^im as a function of 9 for parallel po- 
larization. The forward and backward scattering is recorded 
separately as indicated by the dashed line. Thus the absolute 
intensity of both scans cannot be compared exactly. 

FIG. 3. Theoretical spectra of the THG intensity of a wa- 
ter droplet with a — 10.5 /im in arbitrary units, represented 
linearly in the main part and as a polar plot in the inset. 
Parallel polarization is considered. 
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